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ABSTRACT 

Propagating star formation in a collisional ring galaxy Arp 10 is investigated by a complex ap- 
proach, which includes the broad- and narrow-band photometry, long-slit spectroscopy, and scanning 
Fabry-Perot spectroscopy. The ionized gas velocity field obtained with best spatial resolution to date 
indicates a non-isotropic expansion of the outer ring with a maximum velocity 110 km s _1 . Strong 
vertical and non-circular motions are also seen in the vicinity of the inner ring. Our kinematic data 
suggest that Arp 10 has a small inclination i = 22° and high total (luminous plus dark matter) mass of 
about 10 12 Mq within a 50 kpc radius. The internal extinction in HII regions and extinction-corrected 
Ha fluxes are estimated from the emission lines. The abundance of oxygen 12 + log(0/H) in both 
star-forming rings is about 8.6. The analysis of spectral indices provides an estimate on the propa- 
gation velocities of both rings and metallicity of the pre-collision stellar population. A small "knot" 
near the nucleus of Arp 10, which was previously suspected as a possible candidate for collision, is 
now unambiguously identified as the " intruder" . The intruder spectrum suggests that it was a spiral 
galaxy before the collision and its present mass amounts to at least one-fourth of the total mass of 
Arp 10. A high value of the intruder velocity relative to Arp 10 (~ 480 km s _1 ) suggests that the 
intruder is either gravitationally unbound or marginally bound to the main galaxy. We use a simpli- 
fied two-dimensional hydrodynamic modeling of galaxy collisions to test a collisional origin of Arp 10. 
We confirm that the sizes of the inner and outer rings, maximum expansion velocity of the outer 
ring, and radial profile of the gas circular velocity can be reproduced by a near-central collision with 
\Q i the intruder galaxy, which occurred approximately 85 Myr ago. We acknowledge that an apparent 

CO ' crescent-shaped distribution of Ha emission in the outer ring is caused by a star formation threshold 

', in the gas disk of Arp 10. 
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P ■ give us a bright example of star formation induced by a S e distribution of young stellar clusters (jAppletonl 

radially expanding density waves, which are triggered by H22SD - However, the nature of the pre-colhsion Cartwheel 

C§ i a near-axial passage of a companion galaxy through the ] s stl11 controversial: it might have been either a gas-rich 

disk of a target galaxv dLvnds fcToomreilWfi: HigdoT] low surface brightness galaxy (Korchagm ct al. 2001J) or 



disk of a target galaxy (Lynds & T oomrd ll976: Higdon low HUliace un « llmras & aiax -y ^oicnagm ei ai. zumi) oi 
QUI iCharmandaris fc Appletodll99l ~The expanding a re g.^ lar ^ e E^xy (|Vorobyov fc Bizyaev j 12003) . To 



density waves in the gas disk of a target galaxy leave be- ldentlf y the ^P es ° f dl * S alaxie s that can s ^ 
hind stellar populations of progressively older age. This gemtors to the collisional ring galaxies, we should ex- 
allows us to study the sequence of large-scale star forma- * end th « observational and theoretical studies beyond the 
tion events by considering adjacent locations along the oartwneei. 

galactic radius. A sequential triggering of star forma- The available sample of collisional ring galaxies, i.e. 
tion by expanding density waves makes collisional ring galaxies with an apparent ring morphology and nearby 
galaxies a unique laboratory for studying young stellar companions, is relatively small. However, traces of ex- 
populations pandmg rings and associated star formation can be found 
The best 'studied collisional ring galaxy to date is ™ otherwise regular galaxies. These features include a 
the Cartwheel. Star formation in the Cartwheel is shlft between the peaks of radial distributions of OB- 
concentrated in the two rings, w hich reve al the pres- associations and HII regions m M 51 M 101, Mid some 
ence of expanding density waves fl995l « othe * S alaxies < l S ™ov % Kam berp] | 1976Q re - 
lAppleton fc Struck -Marcelllfl996T) . The non-stationary centh / discovered ring m NGC 922 (|Wong et al. || 2006 |), 
and propagating nature of star formation is manif ested and ring-lik e structure m the inner region of M 31 
by the current radial expansion of the rings (feigdon 1991; ( |Block et al. | |20Q6D . This evidence suggests that head- 
lAmram et all HoM iVorobvov fc Biz vacv 2003), radial 011 collisions with satellites are likely to be a more fre- 

quent phenomenon m the local universe than they were 
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galaxy. We have chosen Arp 10, a galaxy with an appar- 
ent ring morphology, because of its large size and con- 
venient location for our observations (in the Northern 
Hemisphere). Unlike the Cartwheel, Arp 10 shows clear 
signs of the old stellar population in the form of ele- 
ments of spiral arms and prominent bulge. This gives us 
an opportunity to study the development of ring struc- 
tures i n regular spiral galaxies. I n contrast to previous 
work s (|Charmandaris et alj 119931: iMarston fc Appletonl 
1995 ; ICharmandaris fc Appletc^ll996t lBransfor d et al.l 
19981) . we considerably improved the spatial resolution 
of kinematical studies of Arp 10. We also search for 
the companion that is presumably responsible for the 
Arp 10 peculiar ring structure. The new set of observa- 
tional data allows us to create a self-consistent picture of 
collision and subsequent star formation in Arp 10. 

2. OBSERVATIONS 

Observations of Arp 10 were conducted with the 6- 
meter telescope at the Special Astrophysical Observatory 
(Russian Academy of Sciences) and multi-mode f ocal re- 
ducer SCORPIO (see lAfanasiev fc Moiseevl [20051 for de- 
tails). The SCORPIO allows us to carry out the broad- 
and narrow-band photometry, long-slit spectroscopy, and 
scanning Fabry-Perot spectroscopy. The detector is a 
2048 x 2048 pixels CCD EEV42-20. 

2.1. Photometry 

The Ha and Johnson-Cousins B- and R-images were 
taken on September 21/22, 2003. The pixel size was 0.36 
arcsec (with binning 2x2). We obtained several dithered 
images in each broadband filter, with a total integration 
time 600 and 70 sec in the B and R-bands, respectively. 
The mean seeing was 1.6 arcsec. The images were cor- 
rected for bias and flat field, cosmic ray hits, and sky 
background. A saturated nucleus in the B- and R-bands 
was corrected using short expositions. The B-band image 
is shown in Figure [1] 

To obtain the Ha image corrected for the continuum, 
we used the 15 A filter centered at 6720 A for the 
emission line and 15 A filter centered at 6765 A for 
the continuum. The images in these filters were inte- 
grated for 1200 and 600 sec for the Ha and continuum, 
respectively. The Ha image of Arp 10 corrected for the 
continuum and superimposed onto the B-band image is 
shown in Figure [2] 

2.2. Long-slit spectroscopy 

We obtained three long-slit spectra of Arp 10 using a 
1.0 arcsec wide slit. The spectra are taken along three 
cuts, which pass through the center of the galaxy and 
are shown in Figure [2 The summary of spectral observa- 
tions including the supplementary information about the 
spectral range and resolution, corresponding integration 
times, and position angles of the cuts are given in Table[TJ 
The first spectral cut coincides with the major axis of the 
outer ring in Arp 10. The second cut is drawn at the po- 
sition angle 118° and covers a larger spectral range with 
a lower spectral resolution. The third cut passes through 
a small "knot" (presumable satellite-intruder) at 5 arc- 
sec from the center of the galaxy. The spectra were re- 
duced by the standard procedure (bias subtraction, flat 
fielding, night sky subtraction, and wavelength calibra- 
tion), and then flux-calibrated with the help of standard 
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Fig. 1. — B-band image of Arp 10. 



TABLE 1 

Summary of SCORPIO long-slit observations 



# 


Date 


Sp. Range 

A 


Sp. Resol. 

A 


Exp. time 
sec 


P.A. 


1 


Sep. 21/22, 2003 


4150-5890 


5 


6000 


138° 


2 


Nov. 02/03 2003 


3700-7400 


8 


4800 


118° 


3 


Aug. 16/17, 2004 


4040-5600 


5 


2700 


50° 
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40 60 



— Aa, (arcsec) 



Fig. 2. — Ha image of Arp 10 (red) superimposed on the B-image 
(blue). Positions of the long-slit cuts are shown by solid lines. The 
cuts are labelled according to the text. 

stars BD+75 325 and BD+25 4655. The dispersion of 
the spectra along the cuts is 0.71 arcsec pix -1 . The 
seeing was 1.4-1.6 arcsec during the observations. The 
one-dimensional spectra of regions within the outer ring, 
inner ring, and center of Arp 10 are shown in Figure [3l 

2.3. Fabry-Perot Spectroscopy 

The Ha velocity field of Arp 10 was obtained via 
7320 seconds of integration with the SCORPIO in the 
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scanning Fabry-Perot Interferometer (FPI) mode with 
32 spectral channels. Each channel corresponds to 
40 km s _1 . The formal accuracy of our radial veloc- 
ity measurements obtained from the Ha line fitting by a 
Gaussian is 10 km s _1 . The typical procedure of reduc- 
tion o f the FPI frames is described in detail by iMoiseevI 
(|2002f ). The spatial sampling of 0.71 arcsec pix and 
velocity resolution of 125 km s _1 are achieved in our 
frames. The data were smoothed to the spatial resolu- 
tion corresponding to a seeing 2.9 arcsec. 

3. STRUCTURE OF ARP 10 FROM THE B AND Ha 
PHOTOMETRY. 

The broad-band images of Arp 10 indicate the pres- 
ence of an exponential disk and prominent bulge. We 
fit the R-band image w ith a model exponen tial disk 
and bulge (see for details ISmirnova et al .1120061 ) and ob- 
tain the following parameters: the exponential disk scale 
length is 11.4 arcsec and exponential bulge scale length 
is 2.8 arcsec. These scales correspond to 7.18 kpc and 
1.77 kpc, respectively, for the adopted distance of 130 
Mpc (Hq = 70 km s _1 Mpc -1 ). One arcsec corresponds 
to 0.63 kpc and the ratio of the bulge-to-disk luminosities 
is 0.4. 

The Ha emission of Arp 10 comes mostly from two 
star forming rings, which are clearly seen in Figure O 
However, the Ha structure is more complicated than is 
expected from a simplified picture, in which a star for- 
mation wave propagates from the center toward the pe- 
riphery. The outer ring is rather asymmetric, with the 
Ha emission coming mostly from its North- Western part 
(top- right in Figure [2|). At the same time, the inner ring 
has a concentric shape. A long arc of star-forming re- 
gions that traces an outer spiral arm joins the outer ring 
in its South-Eastern part (bottom-left in Figure [2]). A 
similar but a fainter arc is seen on the opposite side of 
Arp 10. 

The structure of the outer disk is best seen in the B- 
band image (Figure [I]) . Outside the outer ring, at least 
three spiral arms can be identified, two of them are re- 
vealed in the Ha image. The diameter of the outer Ha 
ring (28 k pc) is twice smaller than D 2 5 = 57 kpc (ac- 
cording to lde Vaucouleurs et al.lfl99ll ). This is in agree- 
ment with recen t multi-band photome try of 15 northern 
ring galaxies by I Romano et al.l (|2007f ) , who have found 
that stellar disks extend outside Ha emitting rings in 
all galaxies of their sample. Figure [1] suggests that the 
pre-collision Arp 10 might have possessed a complicated 
structure typical for spiral galaxies. 

4. KINEMATICS OF ARP 10 

In this section we derive the kinematics of Arp 10 us- 
ing the long-slit and FPI spectroscopy. Figure 0] shows 
the Ha, [OIII], and stellar line-of-sight velocities for our 
long-slit data (top panel) taken along the first cut. The 
velocities in Figure H are not corrected for inclination. 
It can be noticed that at some radii stars rotate faster 
than the ionized gas, which indicates a presence of es- 
sential non-circular motions in the galaxy. The bot- 
tom panel in Figure |4] demonstrates the stellar velocity 
dispersion distribution, which is typical for very large 
spiral galaxies with high amplitudes of rotatio n curves 
(|van der Kruit k. Freeman] 1 1986t lBottemalll993f ). 

The two-dimensional Ha velocity field of Arp 10 was 
analyzed by a tilted-ring technique in the same manner as 



7 




^ 6 






< 


Center 










o 4 






CP 






S 3 






O 

x 2 


ft Inner ring 




1 


(EW i Outer rinq 
— > U . , . 












4000 5000 6000 7000 

X, A 



Fig. 3. — One-dimensional spectra of regions within the outer 
ring, inner ring, and center of Arp 10 from our long-slit spec- 
troscopy (cut 2). 
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Fig. 4. — Line-of-sight velocities (top panel) of Arp 10 obtained 
with the long-slit spectroscopy. The ionized gas velocities are 
shown by diamonds (H/3) and squares ([OIII]). Stellar velocities 
are designated by filled circles. Radial distribution of the stellar 
velocity disper sion is shown in th e bottom panel. 

described bv lSmirnova et ail (|2006f ). A si milar technique 
was a pplied to the Cartwhe el galaxy by (|Amram et al.| 
1998) and to Arp 10 by ICharmandaris fc Appletonl 
(1996). In contrast to our analysis, Charmandaris & Ap- 
pleton used low-resolution HI data. 

To analyze the velocity field of Arp 10 shown in Fig- 
ure we assume that the line-of-sight velocity of the 
ionized gas V Q b s is a combination of the following four 
components: a systematic velocity V sys , azimuthal veloc- 
ity in the plane of the disk V r j >1 radial expansion velocity 
in the plane of the disk V r , and vertical velocity V z : 

V bs = V sys + V<j> cos <fr sin i + V r sin <f> sin i + 

+ V z cosi. (1) 

Here i designates the inclination of the galactic plane and 
4> denotes the azimuthal angle, which is measured coun- 
terclockwise from the major axis. The position angle of 
an arbitrary point in the disk p. a. is connected with 
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Fig. 5. — Two-dimensional velocity field in Arp 10 obtained with 
the scanning FPL analysis. 

angle 4> as tan(p.a. — P. A.) ~ tan(0) cos(z), where P.A. 
is position angle of the major axis (line-of-nodes) . We 
split the galactic disk of Arp 10 into 48 tilted circular an- 
nuli. In the case of a purely circular motion (P.A.=const, 
i=const, V r = 0, and V z = 0), the best fit position angle 
of the major axis and inclination are 172±2° and 22±4°. 
We obtained these values by masking the inner and outer 
rings. 

The mean rotation curve of Arp 10 is derived from our 
two-dimensional Ha velocity field and is shown by the 
points in the upper panel of Figure Next, we use an 
assumption of quasi-circular motions and calculate P.A. 
independently in each elliptical annulus. The points in 
the middle panel of Figure [6] show P.A. at different dis- 
tances from the center of Arp 10. The points in the 
bottom panel represent a sum of the systematic and ver- 
tical components of the line-of-sight velocity. We note 
that the rotation curve continues to rise slowly even at 
large radii (50 arcsec, or 32 kpc). The dashed lines in all 
three panels of Figure [6] show the values that are used 
to construct the model velocity field using equation (TTJ). 
This field is presented in the top panel of Figure [7J The 
model velocity field was subtracted from the observed 
velocity field to obtain the residual velocity field, which 
is shown in the bottom panel of Figure [7J 

The kinematic data suggest that the best-fit inclina- 
tion is 2 2° rather than the photometric i nclination 49° 
used bv ICharmandaris fc Appletonl (|1996| ). Taking into 
consideration the asymmetric shape of the outer ring in 
Arp 10, this lar ge difference in the v alues of inclination is 
not surprising. iReshetnikovi (|1999l ) noticed that Arp 10 
would not fit into a general Tully-Fisher relation, if an 
inclination of i — 49° is adopted. A lower value of incli- 
nation i — 22° that provides better agreement with the 
Tully-Fisher diagram suggests a total mass of Arp 10 to 
be 7 x 10 11 Mq within a radius of 30 kpc. The value of 
the kinematic position angle of the major axis 172° is 
in poor agreement with the value 125° derived from the 
photometric data (|Paturel et al . 1 120031 ) . Assuming that 
the outer spirals in Figure [T] are trailing, we conclude 
that the galaxy rotates counterclockwise and its Eastern 
part (the left half in Figure [2] and Figure [5]) is closest to 
us. 

The residual velocity field of Figure [7] and radial vari- 
ations of P.A. and V z cos i in Figure [6] indicate the pres- 
ence of two main features. 

1) Velocities in the inner ring suggest significant verti- 
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Fig. 6. — Tilted-rings analysis of the velocity field in Arp 10. Top 
panel: Rotation curve with error bars (circles and vertical bars) 
derived from 2-D velocity field and model rotation curve (dashed) 
used to determine residuals in Figure \E\ Middle panel: position 
angles at different radii and adopted position angle (dashed line). 
Bottom panel: the observed (circles and bars) and model (dashed 
curve) sums of systematic and off-plane components of the line-of- 
sight velocity in Arp 10. 
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Fig. 7. — Top panel: the circular velocity field in Arp 10 con- 
structed using equation JTJ and values derived from the tilted- 
rings analysis. Bottom panel: Line-of-sight velocity residuals in 
Arp 10. A contour map of Ha emission is superimposed by thin 
black curves. 

cal motions of ionized gas with an amplitude of the order 
of 30 km s _1 . The tilted-ring analysis performed for i — 
22° and V r —Q also implies significant vertical velocities 
V z in the vicinity of the inner ring (see middle and bot- 
tom panels in Figure [S]) . 
2) The outer ring expands non-uniformly, with the ex- 
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Fig. 8. — In-plane residual velocities (circles and error bars) 
at the outer ring of Arp 10 at different in-plane position angles 
4> (which are counted from the south counterclockwise). The solid 
curve designates the model expansion velocities obtained in numer- 
ical simulations in § [8] at 85 Myr after the collsion. The dashed 
curve is the same but for 105 Myr. 

pansion velocity in the plane of the disk being larger in 
the North- Western part of the ring (top right in Figure[5]) 
than in the South- Western part. Note that higher expan- 
sion velocities in the outer ring correspond to larger dis- 
tances from the center. This implies that different parts 
of the outer ring has started to expand at the same time 
and from the same position in the disk, which favors the 
collisional origin of the outer ring in Arp 10. 

The radial variations of the kinematic P.A. in the mid- 
dle panel of Figure [6] is a distinctive featu re of radial mo- 
tions, according to iMoiseev et al.l (|2004l ) and references 
therein. Note that we assume that all non-circular mo- 
tions in the outer ring occur in the galactic plane and the 
azimuthal component of non-circular velocities is neg- 
ligible. Strictly speaking, this is true only for the re- 
gions that lie at the minor axis (i.e. regions with P.A. 
= 82° and 262°). In all other regions, the vertical and 
azimuthal components may contribute to the residuals. 

The residual velocities in the inter-ring region may 
achieve rather high values. At the same time, the signal- 
to-noise ratio (S/N) is quite small there, which makes 
the interpretation of gas motions in the inter-ring region 
difficult. 

To see how the expansion velocity in the plane of the 
disk V[ varies along the outer ring, we average V r over 
7 pixel (~ 5 arcsec) circular regions, which are regularly 
spaced along the ring. The resulting distribution of V T 
along the outer ring is shown in Figure [8] as a function 
of angle <j>. Here cj> is measured counterclockwise in the 
plane of Arp 10 from its south direction. The fastest ex- 
panding part (« 110 km s _1 ) of the outer ring at <f> ~ 
130°- 150° corresponds to the highest surface bright- 
ness in Ha. The opposite side of the ring at <f> ~ 310 - 
330° expands considerably slower (up to 30 km s" 1 ). 

5. LONG-SLIT SPECTRA OF ARP 10 

Although we obtained three long-slit radial cuts, only 
the second cut (hereafter, cut 2) is well suited for the 
analysis of emission and absorption lines because of its 
wide spectral coverage and best S/N ratio. 

5.1. Emission lines 

The most intensive emission lines are seen in the star- 
forming rings of Arp 10. Cut 2 in Figure [2] (P.A. =118°) 
covers the regions of most intensive emission and ranges 
from 3700 A to 7400 A, which allows us to estimate 
extinction and heavy element abundances along the cut. 



To estimate the extinction, we use the Ha/H/3 line ratio. 
For the correct extinction and abundance analysis, we 
have to take the absorption lines into account. 

In order to obtain the emission line fluxes, we fit the 
observed spectrum around several selected spectral lines 
([OII]AA3727,3729, H/3, [OIII]A4959, and [OIII]A5007) 
with a combination of model continuum, absorption 
lines, and one Gaussian emission line. In the case of 
blended lines (Ha, [NII]A6548, [NII]A6584), three Gaus- 
sian emission lines were used. The free parameters were 
the width and height of emission lines, and the height 
of continuum. The best-fit parameters were used to ob- 
tain the emission line fluxes. These fluxes at different 
distances from the center of Arp 10 are given in Table [H 

The extinction in the V-band (Ay) from HII regions lo- 
cated at d ifferent radii is e stima ted using the parameter- 
ization of ICharlot & Falll ([2000). The results are shown 
in Table [3] The extinction can be determined only in 
star-forming rings: it attains a mean value of 0.32 mag 
in the inner ring and increases toward the outer ring, 
reaching a maximum value of 4.0 mag in the external 
part of the outer ring. This is in general agreement with 
a propagating star formation scenario - massive young 
stars destroy dust and decrease extinction in the wake 
of a propagating wave of star formation. The observed 
emission fluxes in Table [2] were corrected for extinction 
and then were used to evaluate the oxygen abundance 
in Arp 10. The temperature-sensitive line [OIIIJA4363 is 
not seen in o ur spectra. Instead , we fo llow an approach 
described bv iPilvugin fc Thuanl (|2005l ). which uses the 
excitation parameter P based on the [Oil] AA3727, 3729 
to [OIII]AA4959,5007 relation i? 23 : 



12+log(0/iJ) 



i? 23 + 726.1 + 842. 2P + 337.5P 2 



.96 + 82.76P + 43.98P 2 



1.793i? 23 

(2) 

Where P = Rz/{R 2 + R3), R2 = I[OII]\\3727,372q/ Ih/3, 
R3 = I[OI 1 7]AA4959,5007/ '1 Hf3, and i? 2 3 = i?2 + #3- 

The derived abundance of oxygen will be applied in 
the spectral index modeling in § [5] The shape of absorp- 
tion spectral lines depends on abundances of heavy ele- 
ments. In turn, the emission line fluxes (and hence the 
extinction and abundances) depend on the subtracted 
absorption spectrum. To determine abundances and ex- 
tinction, we use the following procedure. We estimate 
extinction/abundances from emission lines and use the 
resulting values to derive the model absorption spectrum. 
The resulting spectrum in turn is used to derive a new 
approximation for the emission line fluxes and, by im- 
plication, for extinction/abundances. This procedure is 
iterated until convergence. 

The resulting values of P and oxygen abundances 12 
+ log(0/H) are given in Table 03 These values agree 
with a mean value of 1 2 + log(0/H) = 8.6 estimated by 
iBransford et~aT1 {1998) for two HII regions in the outer 
ring of Arp 10. However, w hen a method o f evalu at- 
ing oxygen abundances from IBransford et al.l (|1998l ) is 
applied to our emission line fluxes, the obtained oxygen 
abundances exceed the solar value at all radii. This is 
unlikely for Arp 10 due to its presumably low surface 
density nature (see § 18. ip . We note that the radial oxy- 
gen abundance gradient in Arp 10 is very small, with a 
mean oxygen abundance of 8.55 and 8.64 in the inner 
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and outer rings, respectively. This implies a negligible 
radial metallicity gradient in young stars as well. 

One of the products of our analysis is the radial dis- 
tribution of extinction-corrected fluxes in Ha and H/3, 
which are shown in Table [3] The Ha fluxes were cor- 
rected for [NH]-emission by deblending. We note that 
the [Nil] /Ha flux ratio is remarkably constant along the 
radius of Arp 10. Table [3] indicates that [NII]A6584/Ha 
= 0.37 ± 0.01 and ([NII]A6548 + [NII]A6584)/Ha = 0.48 
± 0.05. These values are typical for quiescent regions 
in regular galaxies rat her than for regions where strong 
shock waves are seen (|Kennicutt k, Keritlll983t iMoiseevI 
2000) . It suggests that propagating gas density waves in 
Arp 10 do not generate strong shock waves. 

5.2. Spectral indices 

The equivalent width of absorption lines in galaxies de- 
pends on a contribution of stellar populations of different 
ages. In case of Arp 10, we assume that stellar popula- 
tions consist of old background stars and young stars. 
The latter were born after a passage of a star forma- 
tion wave associated with an expanding gas density wave. 
We apply several population synthesis codes to construct 
model spectra and estimate the relative amount of old 
and young stars at a specific radius in Arp 10. The prop- 
agating wave of star formation in Arp 10 is approximated 
by a sequence of single-burst events. Spectra of single- 
burst populations are functions of age, metallicity, and 
initial mass function (IMF). We assume that the IMF has 
a universal nature and does not depend on a particular 
location in Arp 10. 

In order to compare model and observed spectra, the 
latter should be corrected for internal extinction. The 
values of V-band extinction from HII regions estimated 
in § 15.11 may differ considerably from the interstellar 
extinction in other regions of Arp 10. Fortunately, 
most spectral indices do not depend on extinction, since 
they are defined in relatively narrow b ands of t h e spec - 
trum. We use index defin it ions by IWorthevi ( 19941 ) , 
Worthev fc Ottavianil (|1997l) . iTrager et alj (jl998l ). and 
Gorgas et al.l (|1999j ). Out of 31 indices, 27 can be mea- 
sured in our spectral range. 

To find the indices that can effectively serve as in- 
dicators of age and metallicity, we use three pop- 
ulation synthesis codes: STARBURST99 (the lat - 
est version is based on [Vazquez fc Leitheren 120051) . 
GAL AXEV tiBruzual fe Chariot! 120031). and Worthev's 
code |WOTHrgj|ll99j). The STA.RBURST99 code pro- 
duces high-resolution synthetic spectra with dispersion 
0.3 A pixel^ 1 in our spectral range. It allows us to 
smooth the model spectra down to our resolution of 5 
A and 8 A. The STARBURST99 code allows to obtain 
spectra of a single-burst population and to choos e param- 
eters for a given IMF. We use the Kroupa IMF (jKroupal 
l2001h and Padova set of stellar models. The spectral 
indices were calculated from the obtained high resolu- 
tion synthetic spectra. The GALAXEV code provides 
spectral indices for a set of ages, metallicities, spectral 
resolutions, and ste llar velocity dis persions. We choose 
the Chabrier IMF (|Chabrierl [20031 ). which is similar to 
the Kroupa IMF, and Padova ev olutionary trac ks. The 
third population synthesis code (|Wortheylll994l ) is con- 
sidered here only for illustrative purposes because of its 
incompleteness at ages of the order of 100 Myr. 




Fig. 9. — Eighteen spectral indices derived for a single-burst 
stellar population with metallicity z = 0.008 (assuming that the 
solar metallicity is 0. 02) at different ages in th ree models: solid 
curves - GALAXEV (Bruzual & Chariot 2003), dashed curves - 
STARBU RST99 HGonzalez Deleado et aLlTT§9a) . dotted curves - 
Worthey HWorthey||1994l ). 

Figure [5] shows how the spectral indices obtained with 
STARBURST99 (dashed lines), GALAXEV (solid lines), 
and Worthey's code (dotted lines) change with age of a 
single-burst population. The metallicity is fixed in all 
models and equal to 0.008 (assuming 0.02 for the solar 
metallicity) , which is close to a mean value of metallicity 
of young population derived for Arp 10 in § 15.11 We note 
that the indices derived with the Worthey's code do not 
extend below 400 Myr. STARBURST99 and GALAXEV 
show general agreement, although a systematic offset in 
the values of indices is seen between the two codes. For 
our analysis, we need spectral indices that show a sig- 
nificant monotonous variation with age. As is seen in 
Figure [H not all indices change monotonously between 
Myr and 500 Myr, which is an expected range of ages 
of young stars in collisional ring galaxies. Some indices 
vary negligibly with age. An additional uncertainty is 
introduced by a dependence of spectral indices on metal- 
licity. In principle, the hydrogen line indices based on the 
strength of H/3, H-y, and H5 lines are the best age indica- 
tors. Unfortunately, a presence of emission lines makes it 
impossible to utilize them. Another age-sensitive index 
Z34000 has to be used with caution, because it is defined 
on a 500 A bandpass in the blue part of spectrum and 
thus depends on the poorly known value of internal ex- 
tinction. For instance, an extinction of Ay = 1 mag 
will decrease D4000 by 7% iGorgas eTai] (|1999l ). Three 
iron-sensitive indices (Fe5406, Fe5709, and Fe5782) are 
not good for the analysis, because their bandpasses inter- 
fere with strong night sky lines in our spectra. Finally, 
we accept 18 spectral indices for the further analysis of 
Arp 10. The observed radial distribution of these indices 
is shown in Figure [TOl 
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Fig. 10. — Eighteen spectral indices (the same as in Figure [9) 
observed in Arp 10 at different distances from the center for long- 
slit cut 2. _ _ 

A brief inspection of Figure [10] and comparison of Fig- 
ure [TU] with Figure [5] show that the indices are sensitive 
to the presence of young stars in the inner and outer 
rings. Systematically smaller values of Fe, Mg, Ca, Ti, 
and D4000 indices are seen in the both rings. Since 
Arp 10 is unlikely to have a similar systematic variation 
in metallicity, the observed variation of indices should be 
interpreted as a presence of young stellar population. 



6. 



SPECTRAL INDEX MODELING OF PROPAGATING 
STAR FORMATION IN ARP 10 



In this section we make predictions about metallicities 
and ring propagation velocities in Arp 10 by fitting the 
radial distribution of model spectral indices to that of 
the observed ones. 

6.1. Model description 

We assume for simplicity that the stellar disk of Arp 10 
consists of young and old stellar populations. The old 
stellar population has an age of either 10 Gyr or 5 Gyr. 
The former value is typical for elliptical galaxies with- 
out ongoing star formation and early- type spiral galaxies. 
The latter value is typical for blue low surface brightness 
and for late-type galaxies. These values are chosen to 
represent two limiting ages of old stellar populations in 
galaxies. The old stellar population in our model consists 
of an exponential disk and bulge, the radial scale lengths 
of which (hold = 7.2 kpc and hi = 1-77 kpc, respectively) 
were obtained from photometry in § [3J The radial sur- 
face density distribution of the old stellar population is 
described as 



^0 



(17.67 exp(-r//i 6 ) + exp(-r//i oW )), (3) 

is derived in § [5] to be 300 M Q pc~ 2 and 
the bulge-to-disk ratio of their central surface densities 



where Eg id 



is 17.67. 

The young population emerges after a passage of two 
circular gas density waves, which propagate from the cen- 
ter of Arp 10 toward its periphery and trigger star forma- 
tion. We assume constant propagation velocities of the 
outer and inner waves (representing the outer and in- 
ner rings) and denote their values by V ou t er and Vi nner , 
respectively. We assume that the radial distribution of 
metallicity of both young and old stars can be approx- 
imated by a linear function of radius. In the following 
text, we normalize metallicities (derived from the spec- 
tral index modeling) to the solar value (z = 0.02) and 
denote the resulting relative metallicities as Z. Using the 
oxygen abundance estimates in § 15.11 and assuming that 
the fractional abundances in Arp 10 are similar to those 
at the Sun, we assume the relative metallicity of young 
stars in the center of Arp 10 to be Zq, v = —0.3 dcx. The 
metallicity gradient of young stars (dZ/dr) y is negligible 
and we set its value to zero. 

We find the metallicity of old stars Z^^ld in the cen- 
ter of Arp 10 and metallicity gradient (dZ / dr) id using 
the least-square fit of model spectral indices to the ob- 
served indices in the inner 0.7 kpc and in the outer re- 
gion between 16.3 kpc and 18.8 kpc. This can be done 
because the region outside the outer ring consists of the 
old stellar population only. The center of Arp 10 does 
contain a certain fraction of the young stellar popula- 
tion, but its contribution is expected to be small due to 
the presence of a high surface brightness bulge. First, 
we use all the indices shown in Figure [10l except for 
CN1, CN2, and NaD which were found to worsen the fit. 
This yields Z 0>old = -0.03 ±0.11 dex and (dZ/dr) old = 
-0.10 ± 0.06 dex kpc -1 . If we consider only Ca, Mg, 
and Ti-dominated indices, we obtain Z 0o i d = 0.08 ± 
0.09 dex and (dZ/dr) oM = -0.08±0.12 dex kpc" 1 . When 
the iron-dominated indices (G4300, Fe4383, Fe4531, and 
Fe5270) are considered, we found that Zo,oid = —0.35 ± 
0.06 dex and (dZ/dr) old = -0.08 ± 0.02 dex kpc -1 . The 
values of (dZ/dr) id derived using different sets of indices 
are similar to each other, whereas the values of i?o,oid ob- 
tained from the iron- and a-element dominated indices 
are significantly different. The values of the abundance 
gradients are t ypical for spiral galax i es (from -0.18 to 0.0 , 
according to Zarits kv et al.l Il994t Ivan Zee et all H99l 
iRolleston et~aIll2000D . 

We assume that the radial surface density distribu- 
tion of young stars can be approximated by an expo- 
nential function E y = Eq exp(— r/h y ). The central sur- 
face density of young stars Eg is expressed in term of 
the central surface density of gas Eg = eE g o, where 
E g o = 22 Mq pc -2 . The coefficient e is a free parameter, 
which characterizes the gas to stars conversion efficiency 
and possible variation of the gas surface density. We de- 
note h y to be the radial scale length of young stars. The 
value of h y is set to be equal or larger than the exponen- 
tial scale length of the old stellar disk (7.2 kpc), since 
the gas disk usually has a shallower radial surface den- 
sity profile than the stellar disk, particularly, in Arp 10 
(|Charmandaris &: Appletonlll996T ). 

We calculate the model spectral indices I m as 

Ur) = ^^, (4) 
if the index is defined in units of equivalent width 
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(|Worthevlll99l or as 



I m (r) = -2.5 log 



B 



old 



B v 



C 



(5) 



if the index is measured in magnitudes. The coefficients 
C, A)id,y, and -Boid.y are determined as 



A 



■old,y 



= E old ^(r)L y (T i d , y ,Z i d:y ) x 



: ^old,y C^old.y) -^old,y) 5 

: S old - v (r)L y (T old , y ,Z i diy ) x 

■0.47 o id )y (T o i d y .Z old y ) 



x 10 

C = £ old (r)MT old ,Z old ) + 



+ ^(r)L v (T y ,Z y ). 



(6) 

The superscripts/subscripts "old" and "y" refer to old 
and young stellar populations, respectively. The V- 
band luminosity of a single-burst population with rel- 
ative metallicity ifoid or Z y and age T D i d or T y is denoted 
as Ly. Spectral indices of a single-burst stllar population 
are denoted as I. 

Our model has four free parameters: the gas-to-star 
conversion efficiency e, scale length h y , V ou t er , and 
Vinner- Chi-square values \ 2 are calculated for all se- 
lected indices using their model and observed values 
(along with their uncertainties) at different radii. An 
integrated value of \ 2 indices is found by summing the 
individual \ 2 indices. We determine the free parameters 
of our model by minimizing the integrated x 2 value. We 
used the Levenberg-Marquardt optimization algorithm. 
The wave propagation velocities may be non-isotropic. 
Therefore, a best match between between model and ob- 
served spectral indices was obtained separately for the 
lower- left and upper-right halves of cut 2 in Figure [H 

6.2. Results of spectral modeling 

A comparison of Figures [9] and [10] indicates that syn- 
thetic spectral indices obtained with STARBURST99 of- 
ten fall outside the range covered by the observed indices. 
This effect does not depend on the S /N ratio of the spec- 
tra or uncertainty in the indices. One of the most promi- 
nent examples of this discrepancy is seen for Mg b: no 
agreement between the model and observations can be 
achieved for any parameters of stellar populations. At 
the same time, the model and observed radial distribu- 
tions of the Mg b index are nearly parallel to each other. 
We find that synthetic indices obtained with GALAXEV 
reproduce the range covered by observed indices better 
than STARBURST99. Therefore, we use the GALAXEV 
code to obtain synthetic indices in this paper. 

First, we consider fifteen observed spectral indices 
shown in Figure 1101 The best match between these in- 
dices and model indices is achieved for the following val- 
ues of free parameters: V oute r — I8O+55 5 ^ m s_1 ' Vinner 
= 46l^ km s _1 , h y — 14 kpc. Although the uncertainty 
in the propagation velocity of the outer ring is rather 
large, the value of V ou ter suggest a very fast expansion 
of the outer ring. The radial distributions of the model 
(open diamonds) and observed (filled squares with error 
bars) indices obtained along the upper half of cut 2 are 
shown in Figure fTTl 

Second, we consider the spectral index D4000 alone. 
Although the D4000 index may be affected by an poorly 
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Fig. 11. — Fifteen model (solid curve and open diamonds) and 
observed (filled squares and error-bars) spectral indices at different 
distances from the center of Arp 10. 

known value of internal extinction, mean gas surface den- 
sity (see § El) implies a rat her small mean value of Ay — 
0.2 mag (|Bohlin et allll97ct ). This value of extinction may 
decrease the D4000 index by only 1.5%. The best match 
between the observed and model D4000 indices is ob- 
tained for the following values: V ou t, 



180t^ 5 km s" 



Vinner = km s , and h y — 65 kpc. Finally, we 

consider two more combinations of spectral indices: Ca, 
Mg, Ti, D4000 and Fe, D4000. The resulting free param- 
eters also suggest a large expansion velocity for the outer 
ring and large value of h y . 

The free parameter e takes a best-fit value of approx- 
imately 0.4. This value is a factor of four larger than 
a typical efficiency of gas-to-stars conversion in star- 
forming regions. A larger value of e suggests that ei- 
ther the surface density of gas is larger than it was as- 
sumed and/or the efficiency of star formation is higher. 
A higher star formation efficienc y was ind e ed rep orted 
for the Cartwheel ring galaxy bv lVorobvovi (|2003f i. The 
above spectral modeling assumed that the old stellar 
population had a mean age of 5 Gyr. A larger value 
of 10 Gyr yields only a slightly worse agreement between 
the model and observed indices. 

7. THE INTRUDER GALAXY 

One of the candidates for the intruder galaxy was in- 
vestigated bv ICharmandaris fe Appletonl (|1996h . They 
found that a companion located 60 arcsec North-East of 
the main galaxy is a distant background galaxy. They 
also suspected another candidate for the intruder galaxy 
which can be seen as a small "knot" 5 arcsec to the South- 
West of the nucleus of Arp 10. We consider this possibil- 
ity and make our third long-slit cut through this object 
(hereafter, the companion). The companion has a typical 
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galactic spectrum with absorption Balmer lines. 

The radial velocity of the companion exceeds the sys- 
tematic velocity of Arp 10 by 480 km s _1 . The position- 
velocity diagram made with the help of absorption lines 
in a wavelength band from 5000 A to 5400 A is shown 
in Figure 1121 The line-of-sight velocity distribution was 
obtained via a cross-correlation of the object spectra 
against a twilight-sky template observed at the same 
night. It is clearly seen that the companion has a notice- 
able velocity dispersion 240 km s , which is comparable 
to that of Arp 10. We consider isodenses that outline 
both Arp 10 and companion in Figure H"2l and estimate a 
maximum range of distances Ar and velocities AV^ cov- 
ered by these isodenses for Arp 10 and companion. The 
kinematic masses of Arp 10 and companion are then es- 
timated using the following formula M = AV 2 Ar/G. 
The resulting target-to-companion mass ratio is about 4, 
which corresponds to the total mass of companion about 
2 x 10 11 Mq. Of course, only the inner part of Arp 10 
was included in the analysis. Arp 10 appears to be much 
more extended than the companion, yet the faint parts 
of the latter might be difficult to observe. 

We note that the target-to-companion mass ratio can 
be larger if the companion is a non-virialized stellar sys- 
tem. On the other hand, the companion may lose a 
substantial part of its mass during the collision (about 
1/3 in a model of lLvnds &; Toomrelll976l ). hence the pre- 
collision mass of the companion can be larger. 

Our B-, R-band and Ha images of Arp 10 do not 
resolve the structure of the companion's outer regions. 
Nevertheless, indirect evidence of the companion type 
comes from the position-vel ocity diagram in Figure 1121 
our H a image, and HI map of lCharmandaris fc Appletonl 
(1996). Neither significant amount of HI nor ionized gas 
are seen in the region of the " knot" . At the same time, 
the stellar velocity dispersion is quite large and the com- 
panion appears to be compact. A visual inspection of 
Figure [T^] indicates that the companion has a size of 5 
arcsec or 3 kpc. This suggests that the companion is lo- 
cated behind the galactic plane of Arp 10 and is shielded 
by the disk of the main galaxy. 

We consider several possibilities for the type of the 
companion, which might have been a disk, elliptic, or 
irregular galaxy before the collision. Disk galaxies are 
expected to contain a modest amount gas, which might 
have been stripped off after the collision. Indeed, a gas 
condens ation can be seen in the HI map o f Arp 10 pub- 
lished bv lCharmandaris fc Appletonl (|1996l ) (see their fig- 
ure 4, top left panel and figure 5). A relatively high linc- 
of-sight velocity of this gas clump suggests that it may 
not belong to Arp 10. The location of the clump does 
not coincide exactly with the companion's nucleus, as is 
indeed expected if the companion is not traveling along 
the line of sight. 

A large total mass of the companion and regular shapes 
of its isophotes argue against a supposition that the com- 
panion was an irregular galaxy before the collision. The 
characteristic age of about 5 Gyr for the companion's 
stellar population (see § 17. ip implies that the companion 
was a spiral rather than an elliptical galaxy. 

7.1. Constraints on the companion's age from its 
spectrum 
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Fig. 12. — Position- Velocity diagram for stellar line-of-sight ve- 
locity distribution (grey-scale) and contour map for the central 20 
arcsec in Arp 10. A difference in the radial velocities of Arp 10 and 
companion is approximately 480 km s . 

The spectrum of the companion shows no emission 
lines. This allows us to use the H/3, H7, and H<5 spectral 
indices to estimate the mean age of stellar population in 
the companion. In Figure [13] we show H/3 versus Mg b 
(upper panel) and H/3 versus Mg 2 (lower panel) rela- 
tions obtained from single-burst models with different 
ages and metallicities. The mo del indices were obtained 
using the GALAXEV package (jBruzual fc Charlotj|2003l) 
and adopted to our spectral resolution and velocity dis- 
persion. The dotted lines in FigurefT3lconnect the indices 
with equal metallicities of 0.008, 0.001, and 0.0008. The 
indices with a solar metallicity z = 0.02 are connected 
by the dash-dotted line. On the other hand, the indices 
with equal ages of 0.4 Gyr, 2 Gyr, 3 Gyr, and 5 Gyr 
are connected by the dashed lines. The solid line with 
squares connects the indices with an age of 10 Gyr. 

The observed indices of the companion are shown by 
the filled diamonds with the error bars in each panel of 
Figure [T31 As can be seen, a single-burst model with 
an age of 5 Gyr and metallicity z = 0.008 yields model 
indices that are similar to those of the companion. An 
estimated mean age of 5 Gyr for the companion's stel- 
lar population implies that the companion was a spiral 
rather than an elliptical galaxy. 

Now, we model the H(3 spectral index under an as- 
sumption of continuous star formation. More specifically, 
we assume that the star formation rate declines exponen- 
tially with time as SFR ~ exp(— ■ t), where 1/(3 is the 
characteristic time scale for sta r formation. In elliptical 
galaxies 1/(3 is less than 1 Gyr (Sandage 1986), whereas 
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Fig. 13. — H(3 versus Mg 2 (top) and H/3 versus Mg b (bottom) 
diagrams. Observed spectral indices of the companion are desig- 
nated by filled diamond with error bars. The dotted lines connect 
spectral indices of a single-burst populations with equal mctallici- 
ties of 0.008, 0.001, and 0.0008. The indices for the solar metallicity 
z = 0.02 are connected by the dash-dotted line. Spectral indices 
of the single-burst populations with equal ages of 0.4 Gyr, 2 Gyr, 
3 Gyr, and 5 Gyr are connected by the dashed lines. The solid line 
with squares connects the indices for the age of 10 Gyr. 

1/(3 — > oo in irregular systems (Gallagh er et al.lll984f) . 
The resulting H/3 spectral index is 

I m = f I{t,z)L v (t,z)SFRdt, (7) 
Jo 

where I(t, z) is the value of the H(3 index for a single- 
burst stellar population with age t and metallicity z, V- 
band luminosity Ly is the same as in equation ([5]), and 
T=15 Gyr. For the estimated metallicity of the compan- 
ion z = 0.008, the observed value H/3 = 2.25 A corre- 
sponds to 1//3 ~ 4 Gyrs, wh ich is typical for Sab galaxies 
(|McGaugh fc de BloklH99l . 

The foreground young population of Arp 10 and possi- 
ble (yet undetected) Balmcr line emission from the com- 
panion may change the values of Balmer indices. This 
may alter the obtained values of mean age and charac- 
teristic time scale for star formation in the companion. 
However, this effect is expected to be small due to a 
lack of gas in the companion. Moreover, the emission- 
independent indices Mg b and Mg 2 in Figure H"3l suggest 
a similar mean age of the companion. 

7.2. Synchronized motion of the outer ring and 
companion 

The companion has a very high velocity of 480 km s _1 
relative to Arp 10. A mean rotation velocity of Arp 10 
is about 300 km s^ 1 (see Fig. [6]), which implies that the 
escape velocity for Arp 10 is about 420 km s _1 and the 
companion may not be gravitationally bound. However, 
a near-central collision suggests that Arp 10 may have 
a massive dark matter halo that extends far beyond 50 
arcsec. In this case, the companion would be marginally 
bound to Arp 10. This might also allow for a much ear- 
lier encounter that might have generated the faint ripples 



et al. 

noted bv lCharmandaris fc Appletor] (|1996( ). The numer- 
ical modeling in § [5] indicates that approximately 85 Myr 
has passed since the collision and the current distance 
to the companion from the center of Arp 10 is about 50 
kpc. At the same time, the projected distance from the 
companion to the center of inner ring of Arp 10 is about 
3 kpc. Hence, the direction of the companion motion 
after the collision nearly coincides with the line of sight. 

The numerical modeling in § [3] suggests that the com- 
panion might have changed significantly its direction of 
motion after the encounter. Position angles of the com- 
panion and the brightest part of the outer ring are dif- 
ferent. The rotation curve of Arp 10 indicates that it 
takes approximately 170 Myr for the brightest part of 
the outer ring to complete a circle. This implies that 
the ring has made about half a turn since its emergence 
85 Myr ago. We reconstruct a synchronized motion of 
the companion and outer ring as follows. The compan- 
ion passed through the disk of Arp 10 about 3 kpc to the 
South-East from the nucleus and generated an expanding 
asymmetric density wave in the disk of Arp 10. The po- 
sition of maximum density perturbation was initially in 
the South-Eastern part of Arp 10 but later relocated to 
the North- Western part due to rotation. The companion 
is now seen 3 kpc to the South- West of the nucleus. 

8. NUMERICAL MODELING OF ARP 10 

In this section we perform a simplified two-dimensional 
hydrodynamic modeling of ARP 10. We want to see if 
the collision with a likely intruder galaxy can reproduce 
the observed sizes of the inner and outer rings and ex- 
pansion velocity of the outer ring. We assume that the 
rings of ARP 10 are a manifestation of gaseous (rather 
than stellar) density waves that propagate in the gas 
disk of the galaxy. This conjecture is motivated by the 
fact that most disk galaxies have their gas disks extend- 
ing far beyond stellar disks. The spectral modeling in 
§ 16.21 suggests that the radial scale length of the stellar 
disk in Arp 10 (« 7.2 kpc) is much smaller than that 
of the gas disk. A neutral hydrogen imaging of Arp 10 
bv lCharmandaris fc Appletonl (|1996l ) also shows that the 
gas disk extends at least 2.7 times the mean size of the 
outer ring. Therefore, we neglect density perturbations 
in the stellar component of Arp 10 and assume that the 
stellar disk serves merely as a source of axisymmctric 
gravitational potential. 

8.1. Model description and basic equations 

Our model galaxy consists of a thin self-gravitating 
gas disk, which evolves in the combined gravitational po- 
tential of a spherical dark matter halo, spherical bulge, 
thin stellar disk, and companion galaxy. The gas disk is 
assumed to be isothermal at T = 9000 K. We use the 
thin-disk approximation to describe the motion of gas 
in Arp 10. In this approximation, the radial extent of 
the gas disk is assumed to be much larger than its ver- 
tical height, discarding the need to solve for the vertical 
motion of gas. A near-vertical impact may send ripples 
through the gas disk of a target galaxy. We assume that 
these vertical oscillations have little effect on the appear- 
ance of Arp 10 due to its low inclination. 

The basic equations governing the dynamics of the gas 
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disk in Arp 10 are 



dv 
~dt 



(v ■ V)v = -V$ g - V$ s 



b,h, 



VP 

~eT' 



(8) 



(9) 



Here, E g is the gas surface density, v is the gas velocity in 
the disk plane, P = c 2 E g is the vertically integrated gas 
pressure, and c s is the sound speed. The gravitational 
potentials of the gas disk <J> g , bulge <&b, stellar disk <I> S , 
dark matter halo $h, and companion $j are provided 
below. 

In order to compute $ g , $b, ^s, and <I>h, we have to 
make assumptions about the total masses and radial pro- 
files of gas, stars, and dark matter in Arp 10. The to- 
tal atomic hydrogen mass in Arp 10 is 1.1 x 10 10 Mq 
(|Charm andaris & Appl etonl 119961 . corrected for a dis- 
tance of 130 Mpc). Hence, we set the total gas mass 
(atomic hydrogen plus helium) within a 50 kpc radius to 
be 1.4 x 10 10 Mq. We neglect a possible contribution 
of molecular hydrogen because its mass is highly uncer- 
tain. The radial distribution of gas in the pre-collision 
Arp 10 is difficult to constrain. We assumed that the 
gas distribution in the pre-collision ARP 10 exponen- 
tially declined with radius and tried different values for 
the central gas surface density E g o and exponential scale 
length r g Q. We found that the Toomre Q-parameter is 
larger or equal 2.0 throughout the disk irrespective of 
E g o and r g o. The Q-parameter may decrease below 2.0 
only if r g o is much smaller than the exponential scale 
length of the stellar disk h \d — 7.2 kpc and the gas 
disk in the pre-collision Arp 10 was strongly concen- 
trated to its center. The latter is unlikely because gas 
disks usually extend further in radius than stellar disks 
(see ICharmandaris fc A pplcton 119961 for the case of 
Arp 10) and thus have shallower radial profiles. This im- 
plies that the pre-collision Arp 10 might have depleted 
its gas content and rea ched a tentative state of gravita- 
tional stability (see e.g. lZasov fc Bizvaevlll996h . We note 
that high values of the Q-parameter are typical for gas 
disks in low surface brightness galaxies. For our simu- 
lations we adopt r g0 = 12 kpc and E g0 = 22 Mq pc 2 . 
We find that these values reproduce best both the radial 
and azimuthal distributions of Ha surface brightness in 
Arp 10. The gravitational potential of a thin gas disk is 
then calculated by (Binney & Tremaine, Sect. 2.8) 



$ g (r, = / r'dr' x 



2tt 



E(r',<p')d<p' 



r' 2 + r 2 — 2rr' cost 



(10) 



This sum is calculated using a FFT technique which ap- 
plies the 2D Fourier convolution theorem for polar coor- 
dinates. 

The present-day rotation curve of Arp 10 is used to set 
constrains on the total masses of stellar and dark mat- 
ter components. The Arp 10 rotation curve flattens at a 
radius of approximately 3 kpc and slowly rises at larger 
radii. We use a steady-state equation of motion for the 
azimuthal component of the gas velocity to fit the model 



and measured rotation curves. We assume that the thin 
stellar disk of Arp 10 has an exponentially declining sur- 
face density profile, whereas the spherical bulge of Arp 10 
has an exponentially declining volume density profile. 
The exponential scale lengths of the bulge hb — 1-77 kpc 
and stellar disk h Q \d — 7.2 kpc are known from the R- 
band photometry of §[3] For the dark matter distribution 
we adopt a modified isothermal sphere 



Ph = 



PhO 



(1 + r/r h y 



(11) 



where pho and are the central volume density and 
characteristic scale length of the dark matter halo. The 
central densities of the stellar disk Eg ld and the bulge pbo, 
and the dark matter halo profile are free parameters. 

We find that the stellar disk and bulge (if considered 
alone) produce a declining rotation curve at radii larger 
than 10 kpc and hence cannot account for the shape of 
the rotation curve in Arp 10, irrespective of the total 
stellar mass. This means that Arp 10 must possess a 
massive dark matter halo. The exact mass partition be- 
tween different components is somewhat uncertain. In 
fact, a strongly concentrated dark matter halo can ap- 
proximately reproduce the Arp 10 rotation curve without 
invoking the need for the stellar disk or bulge. How- 
ever, it would contradict the presence of a considerable 
old stellar population in Arp 10 as implied by spectral 
modeling in § 16.21 After experimenting with different 
values of the free parameters, we find that the rotation 
velocity of Arp 10 at large radii (w 320 km s _1 ) is best 
reproduced by a massive halo with M h = 8 x 10 11 M Q , 
r'h = 1.5 kpc, and pho = 0.6 M pc~ 3 . We also in- 
clude a stellar disk and bulge with masses and central 
densities M s = 1.0 x 10 11 M , Eg ld = 300 M G pc~ 2 
and M b = 4x 10 10 Mq, p h0 = 0.3 Mq pc" 3 , respec- 
tively. All masses are calculated within a 50 kpc radius 
of Arp 10. The masses of the stellar disk and bulge are 
in fact the upper limits, because a more massive stellar 
disk and bulge would generate a declining rotation curve 
at radii larger than 15 kpc, which apparently contra- 
dicts the slowly rising rotation curve of Ar p 10. We note 
that the parameters of the ste llar disk (see iBizvaev et al.l 
l2002tlBizvaev k Kaisinll2004h . when taken together with 
high values of the Q-parameter, strongly suggest a low 
surface brightness nature of the pre-collision Arp 10. It is 
not unexpected that interacting (colliding) galaxies may 
have had low surface brightness progenitors. This is be- 
cause interactions usually produce a new specific struc- 
ture, which is much easier to see if the progenitor galaxies 
have been of low surface brightness. Otherwise, the new 
structure may wash out unless the interaction is very 
strong. The detection of a low surface brightness pro- 
genit or of the interacting galaxy Arp 82 (|Hancock et al.l 
2007) adds more credit to this hypothesis. 

The model rotation curve (solid line), radial distribu- 
tion of the Q-parameter (dashed line) in the pre-collision 
Arp 10, and the measured rotation curve (filled circles) 
of the present-day Arp 10 are shown in Figure [T5J The 
model rotation curve shows good agreement with the 
measured curve except for the inner 3 kpc, where the 
measured curve rises steeper than the model one. We as- 
sume that the steep rise is probably caused by material 
that flows into the inner ring from the inter-ring region. 
Such a radial inflow of gas can indeed be seen in the 
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Fig. 14. — Model rotation curve (solid line) that is assumed in 
hydrodynamic simulations, radial distribution of the Toomre Q- 
parameter (dashed line), and observed rotation curve of Arp 10 
(filled circles). 

model velocity field obtained in the next section. Hence, 
the steep rise of the measured rotation curve could be 
a manifestation of angular momentum conservation and 
the pre-collision Arp 10 might have had no peak at 3 kpc 
and shallower rotation curve in the central region. 

Once the parameters of the stellar disk, bulge, and the 
dark matter halo are determined, we can calculate their 
gravitational potentials. The gravitational potential of 
a thin stellar disk in the plane of the disk, the surface 
density of which declines exponentially wit h radius, can 
be evaluated as (|Binnev fc Tremam"e1ll987l ) 

<£> s = -irGZf d r [h{y)Kx{y) - h(y)K (y)} , (12) 

where y = r/(2h id) and I n and K n are the modified 
Bessel functions of the first and second kinds. The ra- 
dial gravity force per unit mass of a spherical bulge is 
evaluated as 

9£b = GM h (r) 
dr r 2 ' 

where Mb(r) is the mass of the bulge contained within 
radius r. The radial gravity force per unit mass of a 
spherical dark matter can be written as 

= inGpborh [r/r h - arctan(r/r h )] (J^j 

The companion is modeled as a point object with mass 
M c = 2 x 10 11 Mq (see §[7]), which moves in the combined 
gravitational potential of the dark matter halo and bulge. 
The equation of motion for the companion is solved us- 
ing the forth-order Runge-Kutta method with adaptive 
step size control. The gravitational potential of the com- 
panion $i is calculate d using a softened po i nt-ma ss po- 
tential as described in lBinnev fc Tremaind (|1987t ). sect. 
2.8. The softening parameter is set to 5 kpc. 

8.2. Code description 

An Eulerian finite-difference code is used to solve equa- 
tions ([H])-© in polar coordinates (r,<j>). The basic al- 
gorithm of the co de is similar to t hat o f the ZEUS 
code presented by IStone fc Normar] (|1992f ) . The oper- 
ator splitting is utilized to advance in time the depen- 
dent variables in two coordinate directions. The advec- 
tion is treated using the consistent transport method 
of Stone and Norman and the van Leer interpolation 
scheme. The timestep is determined according to the 
usual Courant-Friedrichs-Lewy criterion. The numerical 
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Fig. 15. — Model residual velocity field superimposed on the gas 
surface density distribution at t = 85 Myr after the collision. The 
velocity of 150 km s _1 is indicated in the upper-right corner and 
scale bar on the top is in Mq pc — 2 . The dashed line shows the 
radial cut along which the circular velocity of gas is determined. 

grid has 512 x 512 points, which are uniformly spaced 
in the azimuthal direction and logarithmically spaced in 
the radial direction, i.e., the radial grid points are cho- 
sen to lie between the inner boundary r- m — 0.2 kpc and 
outer boundary r out = 50 kpc so that they are uniformly 
spaced in logr. 

8.3. Results of numerical modeling 

We model the collision by releasing the companion at 
50 kpc above the disk plane of Arp 10. The initial veloc- 
ity of the companion is 490 km s _1 . It falls near axially 
and hits Arp 10 at a distance of approximately 2.7 kpc 
from the center. At the time of impact, the companion 
velocity relative to Arp 10 is approximately 800 km s -1 . 
The gravitational perturbation from the companion sets 
off expanding density waves in the gas disk of Arp 10. 
Figure [T5l shows the residual velocity field superimposed 
on the map of the gas surface density in Arp 10 at 
t m 85 Myr after the collision. The residual velocity field 
is obtained by subtracting the circular velocities from 
the proper velocities of gas. Since we were limited by 
the two-dimensional approach, we could not fine-tune 
the numerical model to reproduce the exact visual ap- 
pearance of Arp 10. Nevertheless, we were able to match 
the mean radii of the inner and outer rings, ~ 3 kpc and 
w 14 kpc, respectively. An overall expansion of the outer 
ring is clearly seen. The model radial expansion velocity 
of the outer ring as a function of azimuth at t = 85 Myr 
is shown by the solid line in Figure [8l The maximum ex- 
pansion velocity is approximately 120 km s , which is 
in good agreement with a measured maximum expansion 
velocity of 110 ± 10 km s _1 . 

The flow of gas in the vicinity of the inner ring is com- 
plicated. The gas appears to stream toward the inner 
ring from the inter-ring region of Arp 10. This inward 
streaming motion of gas provides an extra compression to 
the inner ring, the gas surface density of which becomes 
an order of magnitude higher than that of the outer ring. 
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If outwardly propagating gas density waves are axisym- 
metric, then they can be identified by characteristic sig- 
natures in the azimuthal flow of gas in a galactic disk. 
Indeed, as the wave propagates outward, it pools the 
matter from both sides of the wave toward its current lo- 
cation. This should alter the circular velocity of gas V<f, 
on both sides of the wave due to conservation of angular 
momentum. More specifically, the angular velocity of gas 
that lies ahead of the wave should increase, whereas the 
angular velocity of gas that lies behind the wave should 
decrease. Of course, any considerable deviation from ax- 
ial symmetry of the wave would break the conservation 
of angular momentum. Nevertheless, these characteristic 
signatures can be detected in the radial distribution of 
Vtj,, if the latter is constructed along a narrow radial cut 
rather than being azimuthally averaged. Filled circles 
in Figure [TC] show the radial distribution of V$ obtained 
along the upper half of cut 2 in Fig [2j A characteristic 
double-peaked shape is conspicuous. The outer and in- 
ner rings are located at approximately 15 kpc and 3 kpc, 
respectively. It is seen that is growing across the 
outer ring from 11 kpc to 18 kpc, exactly as predicted 
from conservation of angular momentum of gas that is 
pooled toward the ring. The radial profile of has a 
local minimum at 11 kpc, a radius at which the flow 
of gas is reversed from expansion to contraction toward 
the inner ring (see Figure fT5|) . As a result, V<j> starts to 
grow again at smaller radii and reaches a maximum near 
the position of the inner ring. Inside the inner ring, V$ 
quickly drops to a minimum. Similar double-peak ra- 
dial profiles of Vrf, were found along the radial cuts taken 
at ±10° and ±20° away from cut 2. We note that the 
rotation curve in Figure [6] does not have a well-defined 
double-peaked shape, probably due to azimuthal aver- 
aging along circular rings. We compare the measured 
radial profile of V<f> with that obtained in our numeri- 
cal modeling. That is, we take a radial cut through the 
densest region in the outer ring (dashed line in FigurelTS")) 
and plot the resulting radial profile of V$ along this cut 
by the solid line in Figure [HI Although a quantitative 
agreement is difficult to attain in our simplified thin-disk 
model due to neglected vertical dynamics of Arp 10, a 
qualitative agreement is clearly seen and is encouraging. 
We consider a characteristic double-peak radial profile of 
circular velocity to be powerful evidence in favor of the 
collisional origin of Arp 10. 

A good agreement between the modeled and measured 
sizes and maximum expansion velocities of the rings in 
Arp 10 suggests that approximately 85 Myr have passed 
since the collision. Our simulations indicate that the 
velocity of the companion has dropped to approximately 
480 km s -1 at t = 85 Myr after the collision, which 
is also in agreement with the measured velocity of the 
companion. The first resonant ring forms soon after the 
collision so that the age of the outer ring should be equal 
or less than 85 Myr. The spectroscopic analysis in § [6] 
suggests a similar age, 90 Myr. 

The distribution of Ha emission from the outer ring of 
Arp 10 is characterized by a crescent shape (Figured]). 
We test a conclusion of Ch armandaris et al.l (|1993f ) that 
this peculiar Ha morphology is a result of both the off- 
center collision and star formation threshold in the gas 
disk of Arp 10. We model the Ha su rface brigh t ness i n 
Arp 10 using a method described in IVorobvovl ([2003) . 
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Fig. 16. — Comparison of radial distribution o f th e model circular 
velocities obtained along the radial cut in Figure [T5l (solid line) with 
radial distribution of observed circular velocities along the upper 
half of radial cut 2 in Figure l2l (filled circles). 

The star formation r ate is assumed to obey a Schmidt 
law (|Kennicuttlll998f ) in the form 

S SFR = 2.5x lO -4 ^ 5 , (15) 

where the star formation rate Ssfr and gas surface 
density S gas are measured in units of Mq yr -1 kpc -2 
and Mq pc -2 , respectively. We modify the Schmidt 
law (|15[) by assuming that star formation terminates if 
the Toomre Q-parameter exceeds Q cr — 1.6, a critical 
value for stability of gas disks (characterized by a flat ro- 
tation curve)_a^ainstjocal_jion-axisymmetric perturba- 
tions (|Polvachenko et aT1ll997t ). The star formation rate 
is converted to t he Ha luminosity using a calibration of 
iKennicut t (1983) for solar metallicity 

SFR (M Q yr -1 ) = 7.9 x 10 -42 L Ha (erg s -1 ). (16) 

Metal-poor stars produce more ionizing photons than 
metal-rich stars. The star-forming regions in Arp 10 have 
lower than the solar metallicity and equation (| 1 6 j) is ex- 
pected to give a lower limit on the Ha luminosity. The 
extinction in the gas disk of Arp 10 is taken into account 
by estimating A v from the mod el's known gas sur face 
density via the relation given in iBohlin et all (|l97c? ) for 
i? v = 3.1. The model Ha fluxes are the n corrected for 
the i nternal extinction An a = 0.75 Ay (jCardelli et alJ 
[19891 ). 

Figure [T7] shows the gas surface density distribution 
(left) and Ha surface brightness distribution (right) at 
three consecutive times after the collision t = 45 Myr, 
85 Myr, and 105 Myr. The azimuthal distribution of Ha 
surface brightness in the outer ring has a characteristic 
crescent shape. The values of model Ha surface bright- 
ness Euct in the oute r ring agree approximately with the 
values measured by ICharmandaris et alJ (|1993f ) . More 
specifically, the measured variation of £h q is between 
1.1 x 10 -16 erg s -1 cm -2 arcsec -2 for the faintest part 
of the outer ring and 6.7 x 10 -16 erg s -1 cm -2 arcsec -2 
for its brightest part. These values are not corrected for 
the internal extinction. The model predicts values equal 
or less than 2 x 10 -16 erg s -1 cm -2 arcsec -2 but cor- 
rected for the mean internal extinction in the outer ring 
An a w 1.3 mag as described above. 

Comparison between Ha and gas distributions clearly 
indicates that the Ha emission originates from regions 
with the highest gas surface density. These regions are 
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characterized by Q < Q CI — 1.6, while the rest of the 
outer ring has Q > Q cr — 1.6. Significant azimuthal vari- 
ations of the Q-parameter in the outer ring are easiest 
to reproduce by a slightly off-center collision (2.7 kpc), 
while a bull-eye central collision would require for more 
sophisticated physical conditions in the outer ring (e.g. 
systematic variations in gas temperatures along the ring 
perimeter). Thus we confirm that an off-center colli- 
sion and star formation threshold are most probably re- 
sponsible for the lack of Ha emission from the North- 
Eastern part of t he Arp 10 outer ring , a con clusion pre- 
viously made by ICharmandaris et al.1 (|1993l ) . Our mod- 
eling predicts that the inner ring should be significantly 
less prominent in Ha than the outer ring. This seems to 
contradict the observations. We attribute this disagree- 
ment to the adopted gas isothermality in our numerical 
model - high gas densities in the inner ring might have 
invoked substantial cooling, which had lowered the gas 
temperature and associated star formation threshold. 

It is interesting to consider the model-predicted tem- 
poral variations in the Ha surface brightness shown in 
Figure [T7] (right frames) . In the early phase of evolu- 
tion, Arp 10 was less prominent in Ha, mainly due to a 
higher gas density and associated extinction in the ex- 
panding rings. At the present epoch, Arp 10 is likely 
to have reached a peak in the Ha surface brightness. 
In the next 20 Myr, Arp 10 is expected to fade in Ha, 
mainly due to a drop in the gas surface density of the 
outer ring. The radius of the stellar disk in Arp 10 
at the fiB — 25 mag arc sec" 2 isophote is D25/2 



28 kpc, according to RC3 (|de Vaucouleurs et a.1.1 Il99lh . 
The middle-right panel in Figure [T7] indicates that the 
present Ha ring in Arp 10 is located at about half the 



ra dius of the stellar disk . This tendency was also found 
by I Romano et al] (|2007[ ) for a sample of 15 northern ring 
galaxies. Thus we argue that a preferable location of Ha 
rings at half the stellar disk radius is probably caused by 
high extinction in the early expansion phase (which ob- 
scures Ha emission and hinders detection of small rings) 
and by a drop in the gas surface density below the star 
formation threshold in the later expansion phase. We 
caution that in the real systems there may be other 
mechanisms such as an interaction between rings and spi- 
ral density waves, which affect local star formation and 
Ha fluxes. We also note that the observed outer ring 
in Arp 10 has a noticeably larger eccentricity than the 
model outer ring. This is not a geometry effect because 
of low inclination of Arp 10 (i = 23°). We attribute a 
low eccentricity of our model outer ring to several factors. 
These are the two-dimensional nature of our modeling, 
initially axisymmetric distribution of gas in the disk of 
pre-collision Arp 10, and use of a fixed dark matter halo 
that cannot move in responce to the companion. 

8.4. Global star formation rate 

The gas-to-stellar mass ratio in Arp 10 is about 9% 
within a 50 kpc radius, which is typical for spiral galax- 
ies. The global star formation rate (SFR) in Arp 10 
can be estimated either from far-infrared (FIR) or Ha 
fluxes. The Ha flux is affected by the internal extinc- 
tion, which according to § 15.11 may vary considerably 
across the disk of Arp 10. Therefore, we use the FIR 
fluxes taken from the NED database to estimate the FIR 
luminosity of Arp 10. The FIR fluxes of 0.784 Jy (60 
/jm) and 1.767 Jy (100 /im) suggest the FIR luminosity 
3.84 x 10 10 L Q . This value corresponds to the SFR of 
21 Mq yr -1 , if a conv ersion factor between the FI R lumi- 
nosity and SFR from IThronson fc Telescol (|1986t ) is used 
in combination with a lKennicuttJ (I1983T) initial m ass func- 
tion. We note that according to iZasovl (|1995f ) the con- 
version fac tor can be almost tw i ce low er than the value 
derived by IThronson fc Telescol (1986). Hence, we set a 
lower limit on the SFR in Arp 10 to be 10 Mq yr -1 . 
T his value is still a factor of tw o larger than that derived 
bvlCharmandaris et al.l (|1993l ) using the Ha photometry. 
iRomano et al.l ( 20071 ) also found that SFRs estimated 
from FIR fluxes are also larger than those estimated from 
Ha fluxes. We attribute this difference to the internal ex- 
tinction, which lowers Ha-estimated SFRs. The SFR in 
Arp 1 is very similar to that inferred by iMavva et a.1.1 
(|2005f) for the Cartwheel. 

On the other hand, our numerical simulations predict 
the present SFR in Arp 10 to be about 6 Mq yr -1 . If our 
FIR-based value of 10 - 21 Mq yr -1 is correct, it means 
that the star formation efficiency in Arp 10 is twice to 
four t imes more efficient than that derived by iKennicuttl 
(1998) and used by us in equation (Ti"5"|) . This elevated 
efficiency of star formation appears to be a common fea- 
ture of ring galaxies. For instance, the Cartwheel galaxy 
is also characterized by a star formation efficiency three 
times larger than that of is olated and interacting galax ies 
in the Kennicutt's sample (|Vorobvov fe Bizvaevl l2003). 

Our numerical simulations indicate that the SFR in 
Arp 10 was a nearly linear function of time during the 
last 85 Myr. Hence, the present epoch SFR of 10 - 21 
Mq yr -1 suggests that (4 — 9) x 10 8 Mq of young stars 
were born as a result of the passage of the density wave. 



Propagating star formation in Arp 10 



15 



This is less than 1% of the total stellar disk mass. 

At the same time, Arp 10 has converted 3 - 6% of 
the total gas mass into stars. These values have to be 
certainly taken as lower limits for the star formation ef- 
ficiency because only approximately half of the total gas 
content is located in the inner 16 kpc. 

9. LIMITATIONS OF OUR APPROACH 

The spectral and numerical modeling of Arp 10 per- 
formed by us has several limitations and simplified as- 
sumptions, the impact of which on our results we discuss 
below. An assumption of constant propagation velocity 
of the rings made in § [5] contradicts num erical simula- 
tions (see e.g. IVorobvov fc Bizvaevl [2003^. In the early 
evolutionary stage, the rings propagate faster than in the 
late stage. This tendency can be seen in Figure which 
shows the azimuthal profiles of the gas expansion velocity 
in the outer ring at t=85 Myr (solid line) and t=105 Myr 
(dashed line) after the collision. The assumption of con- 
stant propagation velocity introduces additional uncer- 
tainty to the spectral modeling. 

A faint signal and low-quality of long-slit spectra in the 
inter-ring region and behind the outer ring introduces an 
uncertainty in the spectral index modeling. Because of 
the adopted thin-disk approximation, we cannot model 
the vertical motions in the disk of Arp 10 after the col- 
lision, although these motions are well seen in the ve- 
locity fields in the inner parts of the galaxy. Ideally, a 
three-dimensional hydrodynamic model of Arp 10 should 
be combined with a population synthesis model, and 
the resulting model colors and spectra should be com- 
pared with observational data. High-resolution maps of 
gas density and internal extinction of Arp 10 would be 
needed in this case. 

Nevertheless, the observed kinematics, spectral index 
modeling, and hydrodynamic modeling of Arp 10 are in 
general agreement. This suggests that we envision cor- 
rectly a general picture of the collision and subsequent 
star formation in Arp 10. We postpone a more detailed 
numerical and spectral modeling to a future paper. 

10. CONCLUSIONS 

We apply kinematic, spectral, and numerical ap- 
proaches to develop a self-consistent model of propagat- 
ing star formation in a collisional ring galaxy Arp 10. 
Our kinematic data have the best spatial resolution ob- 
tained for Arp 10 to date. 

We found that the pre-collision Arp 10 was a large 
(Sb?) disk galaxy with an extended low surface density 
exponential stellar disk (scale length 7.2 kpc), substan- 
tial exponential bulge, and massive dark matter halo. 
The spiral arms that are traceable in our B-band and 
Ha images beyond the outer ring suggest that Arp 10 
has had a spiral structure before the collision. Our spec- 
tral index modeling suggests the solar (using all selected 
indices or Mg, Ca, and Ti-related indices) or slightly sub- 
solar (using iron-dominated spectral indices) metallicity 
of the old stellar population in the center of Arp 10. The 
gradient of relative metallicity of the old population is 
-0.08 dex kpc" 1 . 

The available measurements of the total hydrogen mass 



and the shape of the rotation curve suggest high values 
of the Toomre Q-parameter (> 2.6) in the pre-collision 
gas disk of Arp 10, which implies that star formation was 
moderated or even totally suppressed before the collision. 
The Q-parameter may decrease below a tentative stabil- 
ity limit of 1.6 if either a substantial undetected amount 
of molecular gas is present or the gas disk is highly con- 
centrated to the center. Both possibilities appear to be 
unlikely. From the spectra of Arp 10 we estimate the 
oxygen abundance to be 12+log(0/H)=8.6 with little 
variation along the radius. 

Some 85 Myr ago an early-type spiral galaxy (Sab?) 
passed slightly off-center and near-axially through the 
disk of Arp 10. The velocity of the intruder galaxy was 
about 800 km s _1 at the time of closest encounter and 
its mass was about 1/4 the mass of the target galaxy. 
The intruder is now seen as a small "knot" 5 arcsec to 
the South- West from the nucleus of Arp 10. The grav- 
itational perturbation from the intruder generated two 
outwardly expanding density waves, which are presently 
seen in the B-band and Ha images as bright asymmetric 
rings. The expansion velocity of gas in the outer ring is 
a strong function of azimuth, with a maximum value of 
about 110 ± 10 km s _1 in the brightest part of the ring 
and a minimum value of about 25 ± 5 km s . The inner 
ring has a more symmetric shape than the outer ring and 
expands with a velocity less than 50 km s _1 . The prop- 
agating density waves in Arp 10 do not generate strong 
shock waves. 

Star formation is triggered at the crest of an expanding 
gas density wave, if the gas surface density exceeds a star 
formation threshold. An apparent crescent-like distribu- 
tion of Ha emission in the outer ring is most probably 
caused by a variation of the Toomre Q-parameter along 
the ring. The current star formation rate is estimated 
from the far-infrared flux to be 10 - 21 M Q yr" 1 . The 
expanding waves of star formation have created less than 
1% of the total stellar disk mass since the collision. 
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TABLE 2 
Emission lines in Arp 10 



r, kpc 



H/3 



Fluxes in lines, in units of 10 
Ha [Nil] 6548 [NH]6584 



erg s cm 
[OH] 3728 [QUI] 5007 



[OIII]4959 



-1.36 
1.36 
11.59 
14.54 
16.36 



8.52±0.34 
8.06±0.35 
4.60±0.10 
6.15±0.09 
2.67±0.05 



27.6±0.32 
23.1±0.28 
11.3±0.14 
21.4±0.23 
17.2±0.19 



2.93±0.12 
4.43±0.18 
1.10±0.04 
1.90±0.03 
1.37±0.04 



9.66±0.15 
8.63±0.14 
4.10±0.07 
7.98±0.10 
6.70±0.09 



2.74±0.23 
2.36±0.26 
4.53±0.19 
4.48±0.15 
1.19±0.09 



0.48±0.12 
0.24±0.06 
0.85±0.04 
2.51±0.05 
1.23±0.04 



0.22±0.03 
0.34±0.02 
0.23±0.02 



Note. — Emission line fluxes at different radii r in Arp 10. 



TABLE 3 

Extinction and abundances in Arp 10 



r 

kpc 


P 


A v 

mag 


12+log(0/H) 
dcx 


fcor(Ha) 

10 -15 erg s _1 cm -2 


fcor(H/3) 

erg s _1 cm -2 


-1.36 


0.132±0.018 


0.64±0.03 


8.57±0.22 


4.63±0.34 


1.61±0.34 


1.36 


0.091±0.017 


O.OOiO.OO 


8.53±0.28 


2.31±0.28 


0.81±0.35 


11.59 


0.190±0.011 


O.OOiO.OO 


8.53±0.10 


1.13±0.14 


0.46±0.10 


14.54 


0.336±0.010 


1.01±0.02 


8.65±0.06 


4.88±0.25 


1.70±0.09 


16.36 


0.320±0.018 


4.25±0.09 


8.62±0.10 


54.6±1.16 


19.0±0.39 



Note. — Extinction Ay in HII regions, excitation parameter P (see text), oxygen abundance 12+log(0/H), and Ha and H/3 fluxes 
corrected for the extinction (f cor (Ha), f cor (H/3)) in units of 10~ 15 erg s _1 cm~ 2 at different radii r in Arp 10. 



